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The integrated Sachs- Wolfe (ISW) effect in the cosmic microwave background (CMB) as measured 
through its correlation with galaxies provides a unique opportunity to study the dynamics of the 
dark energy through its large scale clustering properties. Ultimately, a deep all-sky galaxy survey 
out to 2 ~ 2 can make a ~ lOcr or ~ 10% measurement of the correlation and limit ~ 3% changes 
in the gravitational potential or total density fluctuation due to dark energy clustering on the Gpc 
scale. A canonical single scalar field or quintessence model predicts that these clustering effects will 
appear on the horizon scale with a strength that reflects the evolution of the dark energy density. 
In terms of a constant equation of state, this would allow tests of the quintessence prediction for 
models where |1 -I- w| > 0.05. 



I. INTRODUCTION 

The two generic signatures of a dynamical form of 
the dark energy are an energy density that evolves as 
a function of time and one that varies as a function of 
space. Clustering in the dark energy is in fact an un- 
avoidable consequence of time evolution on superhorizon 
scales. Prospects for measuring the density evolution or 
equation of state of the dark energy are well known. De- 
tecting its spatial clustering is considerably harder but 
of at least equal fundamental importance. In particular, 
the clustering tests the quintessence hypothesis, dark en- 
ergy as a single scalar field with a canonical kinetic term 

m 

Perhaps the best hope of studying the clustering of 
the dark energy lies in its effect on the cosmic microwave 
background (CMB). In general, the effect of having a 
"smooth" component that alters the expansion rate with- 
out clustering with the dark matter is to make gravita- 
tional potentials decay in the linear regime. As a CMB 
photon transits a decaying potential well, it picks up a 
net blueshift. When integrated across the potential hills 
and wells along the line of sight, this effect leaves a tem- 
perature anisotropy across the line of sight called the in- 
tegrated Sachs Wolfe (ISW) eff'ect [IH- On scales where 
the dark energy is clustered, this effect is reduced in a po- 
tentially measurable way. It is a curious coincidence that 
the observed large angle temperature anisotropy is also 
suppressed com par ed with a cosmological constant model 
(e.g. [3 and [gI l2lj for a possible ISW connection). 

Since the decay of the gravitational potential is a direct 
consequence of the dark energy, the ISW effect is more 
sensitive to changes in the clustering of the dark energy 
than other clustering statistics. Unfortunately this sen- 
sitivity is hidden in the CMB by the anisotropy from re- 
combination. The ISW contributions can be isolated by 
cross correlating the temperature field with other tracers 



of the gravitational potential, in particular galaxies 
and gravitational lensing statistics |^ 0, 0| . 

Recently the ISW effect has been detected by cross 
correlating temperature maps from WMAP with radio 
galaxies, optical galaxies, infra-red galaxies and X-ray 
sources ^2 Q| ■ The current measurements are 

not sufficient for much beyond a simple detection of the 
effect, but do indicate that the technical issues involved 
are well within in the capabilities of current and future 
surveys. 

In this paper we study the prospects for constraining 
the smoothness of the dark energy with future galaxy 
surveys. This study extends that of Afshordi TS] who 
mainly considered prospects for constraining the equa- 
tion of state of the dark energy under the small scale 
Limber approximation. ISW constraints on a smooth 
equation of state are not competitive with other probes 
of dark energy evolution due to the unavoidable cosmic 
variance of the CMB temperature field, whereas they 
may provide the best means of studying the clustering 
of the dark energy. 

The outline of the paper is as follows. In we re- 
view the predictions and parameterization of the large 
scale clustering of the dark energy. In ijllll we calcu- 
late the ISW-galaxy correlation under the general all-sky 
framework applicable to correlated galaxy distributions 
and fine redshift resolution. In ^IIVI we study the poten- 
tial of galaxy surveys to constrain the clustering of the 
dark energy. 



II. DARK ENERGY CLUSTERING 

With adiabatic initial conditions and in the absence of 
non-gravitational forces, all energy density components, 
including the dark energy, are clustered. In this limit, the 
gravitational potential or Newtonian curvature $ evolves 




FIG. 1; Fractional difference in the gravitational potential 
today between smooth and clustered dark energy models as a 
function of w (upper panel) where Ode and h are adjusted to 
keep the distance to recombination and the expansion rate at 
high-2 fixed. Also shown are the predictions for quintessence 
dark energy (canonical kinetic term with sound speed Ce = 1) 
near the horizon scale k — 10"* Mpc~^ and near the broad 
peak of the ISW effect at low multipoles k = 10~^ Mpc~^. 



in a flat universe according to the relation [l^ ll^ 
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where Q is the initial comoving curvature and H is 
the Hubble parameter. Note that during epochs when 
the expansion is dominated by a species with a con- 
stant equation of state parameter wt = Pt/ Pt, H{a) cx 
Q,-3(i+u>r)/2 ^Yxe gravitational potential is constant 
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With only gravitational forces, the ISW effect vanishes 
in a flat, constant wt universe except for the special case 
of Wt = — 1 where the constant is zero. 

To accelerate the expansion, dark energy requires rel- 
ativistic stresses where the pressure is comparable to the 
energy density in magnitude. Stress gradients can pre- 
vent the clustering of the dark energy on small scales. 
Given dark energy that is smooth compared with the 
dark matter, the gravitational potential evolves in a flat 



universe as 
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FIG. 2: Fractional change in the gravitational potential in 
the transition regime between smooth ($3) and clustered dark 
energy for several choices of the sound speed and & w — —0.8 
model. The quintessence prediction corresponds to = 1. 



Here w{a) = pde/pde is the equation of state for the 
dark energy and f7DE(a) = PDE(a)/Pcrit(a)- Where no 
argument is given a = 1 is assumed and the Hubble con- 
stant Ho = H{a = 1) = 100/ikm s~i Mpc"^ The initial 
conditions for scales well above the horizon at matter- 
radiation equality are provided by Eq. JQ) with a starting 
epoch in the matter-dominated limit. For smaller scales, 
the initial conditions are modified by the usual transfer 
function to account for radiation stresses. 

In Fig. n] (upper panel), wc show the fractional differ- 
ence in the gravitational potential today for the smooth 
and clustered regime relative to a cosmological constant 
w — —I model. In this model $ has decayed from 
— in the matter dominated wt = regime 

to <I>A = 0.75<i>m by the present. 

On the right axis, we normalize the curves relative to 
the net change in the gravitational potential $,„ — <I>a = 
0.25<I',„ in the w = —I model. Since the ISW effect 
is sensitive to the change in the gravitational potential, 
this factor of 3 enhancement of the difference reflects the 
observational effect. This enhancement is the reason why 
the ISW effect is more sensitive to dark energy clustering 
than other measures of clustering. 

Note that here and throughout when considering mod- 
els of different w, Que and h are adjusted to keep the 
well-constrained distance to recombination and high red- 
shift expansion rate fixed (Fig ^ lower panel). This 
adjustment guarantees that the CMB predictions for 
the various dark energy models will be indistinguishable 
aside from the ISW effect of interest (see also Fig. O. 

The transition between the smooth and clustered 
regimes depends on the physical model for the dark en- 
ergy. It is usefully parameterized by the effective sound 
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speed of the dark energy Ce, defined through 
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in a "rest frame" coordinate system where the momen- 
tum density of the dark energy vanishes In terms of 
the gravitational potential phenomenology, the minimum 
of the sound horizon and (particle) horizon separates the 
smooth and clustered regime. Above the horizon, there 
is no unique definition of smoothness for a dark energy 
component whose density varies with time given the free- 
dom to choose the time slicing in general relativity . 

In a scalar field dark energy model, the sound speed 
so defined is exactly given in linear theory by the form of 
the kinetic energy as a function of the field [isl 0| . For 
a canonical kinetic term, the sound speed Ce = 1. As a 
dark energy candidate, a canonical scalar field is dubbed 
quintessence and one with an alternative kinetic term 
k-essence |20|. For quintessence, the transition between 
the smooth and clustered regimes thus occurs near the 
horizon scale. In Fig. ^ we also show the gravitational 
potential amplitude for a scalar field model at k = 10~^ 
Mpc~^ and k = 10~^ Mpc~^. Note that the horizon scale 



'70 



^ dt 



14Gpc , 



(5) 



and so the horizon wavenumber kn = i^q^ ~ 7 x 10~^ 
Mpc~^. In the w = —1 model contributions to the 
ISW effect at the scale of the quadrupole are broadly 
distributed around k = 10^^ Mpc~^ due to projection 
effects (see e.g. [23 )■ Thus a fair fraction of the contri- 
butions to the ISW effect will come from the transition 
regime if u; ^ -1 (e.g. HIISIll). 

In Fig. 121 we show the effect of dark energy clustering 
as a function of scale for Cg — 10^^'^, 10^^ . . . 10^ for 
w = —0.8 the smooth limit of Eq. ||2Jl. We employ 
these numerical results for the potential spectrum at the 
present in the following sections. For the time evolution 
employed in the next section, we use the interpolation 
between solutions of the clustered and smooth regimes 
in Eq. (P) and (jJl given in j^, vahd for Ce < 1. 

Note that the clustering effects saturate close to Cg ~ 1 . 
For this reason, we will take models with Cg < 1. For 
illustrative purposes we will also typically employ the 
w = —0.8 model but discuss the conversion of constraints 
to general dark energy models in 3IVI 



III. ISW CROSS-CORRELATION 

The effects of dark energy clustering can be seen on 
the two-point correlation between observed fields that 
depend on the gravitational potential such as the CMB 
temperature and the galaxy number density. Let us as- 
sume that a field x{n) as a function of the angular posi- 
tion n on the sky is a weighted projection of the potential 
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where the weight W'" can include differential operators 
acting on the field. Here D = dz/H{a) is the co- 
moving distance to redshift z. Then the angular cross 
correlation between two observed fields x and x' is given 
by 

(x(n)x'(n')) = P'i^ ■ ft')^C^r' , (7) 



where Pp is the Legendre polynomial and C^^ is the cross 
power spectrum. The angular power spectrum itself is 
given by 

/d'^k , 
^I^{k)I^ {k)P'''^ik;0), (8) 

where P**(fc; z) is the 3D power spectrum of the poten- 
tial field 

($(k; 0)$(k'; 0)) = (27r)3,S(k - k')P**(fc; 0) (9) 
and the weights 
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Note that in spite of the fact that in reasonable dark 
energy models the gravitational potential changes most 
rapidly around the current epoch, the local contributions 
to the angular power spectrum are strongly suppressed 
by the ji projection factors. The contributions from fi- 
nite k locally appear in the monopole moment of the 
temperature field which makes a small and inseparable 
contribution to the background temperature. 

In Fig. Owe show CMB temperature power spectra for 
the w = -0.8 (f^DE = 0.68) model with the ISW contri- 
bution separated out for c,, — 1 and Ce = 0.1. In addition 
to the dark energy parameters, we will assume a physical 
non-relativistic matter density of 0,mh? = 0.14, a physi- 
cal baryon density of Q.i,h^ — 0.024, a reionization opti- 
cal depth of r = 0.17, an initial scale invariant comoving 
curvature spectrum with amplitude 5q = 5.07 x 10~^ 
(WMAP A = 0.87 [23 ) and slope n = l. 

Although the ISW contributions for the two models in 
Fig. 121 differ by up to a factor of 2 as is consistent with 
Fig. n the difference in the total temperature field is 
considerably less due to contributions from other effects. 
Given the cosmic variance of the temperature field also 
shown, the subtle difference is difficult to detect from the 
CMB alone. 



4 



10-9 



10-11 



10-12 











r 






total ] 






N \ 
N \ 
N \ 
\\ 


Cg=l 

c,=0.1 


w=-0.8 


\\ 
\\ 








10 

/ 


100 



FIG. 3: CMB temperature power spectrum for the lo = —0.8 
model compared with the WMAP 1 year data with noise er- 
rors and the 68% and 95% cosmic variance confidence bands 
plotted for Ce = 1. The Ce = 0.1 model is difficult to distin- 
guish from the Ce = 1 model even with perfect data despite 
the fact that the ISW effect changes by up to a factor of 2. 
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FIG. 4: Bias and redshift distribution of the galaxies. The 
total number density n^*°' (here multiplied by 2/n^'°', with 
TiA ~ 70 gal arcmin^^ for clarity) defines the linear bias 
(0; z) under the halo model [23 |. This total number den- 
sity is divided into photometric redshift bins of 5cr(z) with 
photometric redshift errors of ai^z) = 0.03(1 + z). 



The intrinsic sensitivity of the ISW effect to dark en- 
ergy clustering can be extracted from cross correlations 
with the galaxy number density fluctuations. For the 
galaxy clustering, there is a number density fluctuation 
field for each population of galaxies. We denote the zth 
population as having a number density fluctuation gi. 
The weights are given by 

W^^^fc; z) = ^ ^ f AV « fo.. (fc. ,) , (12) 

where the n^* (z) is the redshift distribution of the popula- 
tion with an angular number density of J dz rfi^ (z) = n\ 
and W^(k\z) is the galaxy bias of the population. For 
the galaxy bias we take the parameterized halo model 
described in . As long as galaxy number fluctuations 
trace the total density fluctuation or gravitational po- 
tential in the linear regime, the details of this model do 
not affect the end results. This is because in the lin- 
ear regime, the galaxy auto correlation C^'^^ determines 
the bias parameters for the ISW-galaxy cross correlation. 
Marginalization of halo parameters mainly eliminates the 
small amount of information coming from the non-linear 
regime. 

For illustrative purposes, let us take populations to be 
selected from an overall redshift distribution of the form 



(z) cx z 



(13) 



with Zn adjusted to give a median redshift of z = 1.5 and 
a total galaxy number density of n^* = 70 gal/arcmin^. 
Under the halo model of 23] this number density deter- 
mines the linear bias shown in Fig.^ These specifications 



are close to what can be achieved from an LSST type sur- 
vey ^3 . This parent distribution can be subdivided into 
multiple populations 



'(z) = ^n^'^(z) 
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through photometric redshifts. To approximate the red- 
shift binning, let us suppose that photometric redshift 
estimates are distributed as a Gaussian with an rms fluc- 
tuation of cr(z). A top hat cut in photometric redshift 
then becomes smooth overlapping distributions in actual 
redshift 



(z) = in^'°'(z) 
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where erfc is the complementary error function. We take 
as a fiducial model a{z) — 0.03(1 + z) and choose the 
binning to span 5(t(z) for a total of 5 bins out to z = 1 
and 10 bins out to z = 3 (see Fig. 0)). 

In Fig. 13 we show the cross power spectra of these 
binned galaxy fields with the ISW effect in the upper 
panel for w = —0.8 and Cg = 1 and in the lower panel the 
ratio of spectra for Ce = 0.1 and this model. Note that 
due to projection effects, the correlation in the higher 
redshift bins peak at higher (. or smaller angular scales 
where the cosmic variance is smaller. 

The amplitude of the effect on the cross correlation as 
parameterized by Ce depends strongly on the background 
equation of state. In Fig. |H| we show the cross correla- 
tion of the z = 0.75 — 1 galaxy bin for w = —1, —0.8, 
—0.6, —0.4. Note that the variations between the models 
at Cg = 1 are relatively small (upper panel) whereas the 
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FIG. 5: Galaxy-CMB cross correlation through the ISW ef- 
fect for the w = —0.8 model. Upper panel: Ce — 1 cross power 
spectra for 10 different galaxy redshift populations of Fig. 2] 
Lower panel: ratio of the Ce = 0.1 and Ce = 1 power spectra. 
Note that the intrinsic 10% change in the gravitational po- 
tential has been amplified by up to a factor of 3 through the 
sensitivity of the ISW effect to changes in the gravitational 
potential. 



difference between Ce = 1 and Ce = 0.1 increases as w 
increases. The consequence is that the cross correlation 
is a relatively insensitive probe of the dark energy equa- 
tion of state given the fixed angular diameter distance 
to recombination and is mainly useful as a probe of dark 
energy clustering. 



IV. CLUSTERING FORECASTS 

Interpreting forecasts for the detectability of dark en- 
ergy clustering requires special care for several reasons. 
Aside from the special case of a cosmological constant, 
no physical model can keep the dark energy smooth near 
or above the horizon scale. It is therefore not physically 
meaningful to compare a given dark energy model with 
and without dark energy perturbations. Although the 
sound speed parameterizes the physical scale of the tran- 
sition to smoothness below the horizon, the amplitude of 
the transition depends strongly on the dark energy equa- 
tion of state w for adiabatic initial conditions or track- 
ing models which lose their dependence on initial fluc- 
tuations. Since changing the sound speed can shift the 
effect outside of an observed range, differences in mod- 
els around a point in the parameter space do not reveal 



FIG. 6: Galaxy-CMB cross correlation as a function of w with 
the high-2 universe fixed as in Fig. The cross correlation 
itself is fairly insensitive to the equation of state for fixed 
sound speed (ce = 1 upper panel) but becomes increasingly 
more sensitive to the sound speed as w increases. 



global constraints. 

For these reasons, constraints involving the sound 
speed should be translated into statements as to how 
smooth the dark energy is out to a given physical scale. 
Such statements will also retain greater validity outside 
of the adiabatic and constant sound speed class of mod- 
els. 

With this in mind let us begin with the usual Fisher 
approach to parameter estimation forecasts. Consider 
the a survey which makes noisy measurements of a set of 
fields Xi with a noise power spectra of N^^^^ . Then the 
total, signal plus noise, power spectra are given by 
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(16) 



Given a parameterization of the signal power spectrum 
with a set of parameters Pa , the information in the survey 
on these parameters is quantified by the Fisher matrix 
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where the sum is over bands of width Al in the power 
spectra and /sky is the amount of sky covered by the 
survey. Here we have suppressed the (xi , Xj ) indices in a 
matrix notation and 
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FIG. 7; Total significance of tfie cross correlation detection 
for galaxies between z = and Zmax and multipoles ^min < 
£ < 1000. Half the detection significance is supplied by z < 1 
galaxies and most of the information comes from 10 < f < 
100. 



where the derivatives are evaluated at a fiducial model 
tlirough a finite difference approximation. The inverse 
Fisher matrix approximates the covariance matrix of the 
parameters w (F^^). We will take as a fiducial model 
the w = -0.8 model of Fig. El 

For the CMB temperature field 8, E-polarization field 
E and galaxy fields in question, we will assume white 
noise power spectra 
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where UAi is the galaxy number density in bin i. All noise 
cross power spectra are assumed to vanish. For the CMB, 
we will take for illustrative purposes noise specifications 
close to the Planck satellite Ar = Ap/-\/2 = 40/iK- 
arcmin and 6'fwhm = 7 arcmin. The details of this noise 
choice are not critical as long as the temperature field is 
sample variance limited out to £ ~ few xlO^. For the 
galaxies, we use the specifications given in the previous 
section but again the number densities are sufficiently 
high so as to make the shot noise subdominant for the £ 
range in question. 

The total signal to noise in the cross power spectra can 
be quantified through a simple one parameter family for 
the spectra 
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such that the Fisher matrix in A gives the significance of 
the detection given a fixed template for the shapes and 
relative amplitudes of all the cross power spectra 



total 



(21) 



Note that here the noise includes the sample variance of 
the fields as it represents the signal-to-noise ratio of the 
correlation detection. 

This total significance is shown in Fig. |7| and is ap- 
proximately 100, or a signal-to-noise ratio of 10 for the 
full range of Zmax < 3 and £ = 2 — 1000 for a full sky 
survey /sky = 1. Even with multiple power spectra from 
galaxies across the whole range of the acceleration period, 
the signal-to-noise remains moderate due to the sample 
variance of the single CMB temperature field TE^. 

Note that for the w = —0.8 fiducial model roughly 
half of the information comes from galaxies at z > 1 
which are accessible only to a very deep survey. (This 
contribution decreasess somewhat as w — > — 1 leading to 
a slightly smaller total signal to noise ^3-) On the other 
hand, doubling the information only increases the signal 
to noise by V2- Furthermore, most of the information 
comes from intermediate multipoles 10 < £ < 100 which 
are accessible to a deep but not necessarily all-sky survey. 
With increasing redshift, the weight moves to higher £ at 
higher z due to projection effects. 

The total signal-to-noise ratio gives a rough quantifica- 
tion of prospects for dark energy clustering constraints. 
A S/N of 10 implies that the amplitude of the cross power 
spectrum can be measured to about 10% from the inter- 
mediate multipoles. Since the ISW effect is sensitive to 
the change in the potential, it gains a factor of ~ 3 en- 
hancement in sensitivity to potential variations. Thus 
roughly the cross correlation can detect a 10% variation 
in $ as in the w = —0.8 model at ^ 3a. This £ range 
implies that this ability applies to a transition scale of 
1% - 10% of the horizon. 

This order of magnitude estimate is borne out by a 
more quantitative treatment. Here we employ the 7 cos- 
mological parameters that determine the fiducial model 
in Fig. 131 as Fisher parameters. To account for the dark 
energy clustering we add the sound speed as a parameter. 
Because the effect of the sound speed is only apparent for 
variations that span an order of magnitude in Cg, we take 
this parameter to be 



log 
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where A logj^Q Ce is the step away from Ce — 1 used to 
define the parameter derivative in Eq. (|18|l . This nor- 
malization factor reflects the fact that the Fisher errors 
should be interpreted as the significance of the observed 
difference between the actual models used to define the 
derivative 
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FIG. 8: Significance of tfie separation between a given Ce 
model and a quintessence model of Ce = 1 given an all sky 
survey of galaxies and allowing for parameter degeneracies 
around a w = —0.8 model. Note that in comparison with 
the total detection significance, dark energy clustering model 
separation requires lower multipoles and utilizes more high-z 
galaxies as Ce ^ 1 and the effects are confined to the largest 
scales. 
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FIG. 9: Conversion of the sound speed constraints into more 
robust dark energy smoothness constraints. Shown are the 
projected Icr errors on the change in the gravitational poten- 
tial due to dark energy clustering as a function of scale. The 
galaxy-ISW correlation can ultimately constrain the smooth- 
ness of the dark energy at the IGpc scale to ~ 3% in the 
potential. 



accounting for parameter degeneracies. 

To these cosmological parameters, we supplement 5 ad- 
ditional "nuisance" parameters per galaxy redshift bin 
that define the galaxy bias. These parameters, given ex- 
plicitly in , define the occupation of galaxies in dark 
matter halos. We then have a total of 58 parameters. 
Though marginalization of these extra parameters can in 
general degrade constraints due to parameter degenera- 
cies, in this case the information in the auto power spec- 
tra of the CMB and galaxies suffices to make the con- 
straints quite comparable to those from employing the 
cross power spectrum only but with the other parame- 
ters fixed. The addition of the galaxy auto correlation 
power spectrum actually enhances the significance of the 
separation between Cg = 1 and Cg <C 1 models and begins 
to play a role at Cg ^ 0.1. Similar to massive neutrinos, 
for very low sound speeds, the transition occurs on scales 
where the galaxy power spectra are well measured and 
the ISW effect is negligible. However there may be sub- 
tleties in the Cg 1 regime involved with galaxy bias as 
to whether the galaxies trace the total density perturba- 
tion or that of the dark matter. 

We show the significance of the separation between 
models with different sound speeds in Fig. [S] including 
marginalization over w and other parameters. Recall 
that for the w — —0.8 fiducial model shown here the 
amplitude of the effect is ~ 10% in the gravitational po- 
tential. The model with Cg ~ 0.1 is distinguished at a 
SjN K, 2.5 for /sky = 1. Note however that the signifi- 
cance drops off rapidly as Cg ^ 1 as the changes become 
confined to the low multipoles. Thus one would draw an 



incorrect inference if one employed log^o Cg directly as a 
Fisher parameter and let the derivative be approximated 
by finite difference with A logj^Q Cg 0. 

Utilization of the ISW-galaxy correlation for dark en- 
ergy clustering does shift the optimal redshift and multi- 
pole range to higher redshifts and lower multipoles com- 
pared with a simple detection significance criteria. Be- 
cause the effects are confined to large physical scales, 
they can only be observed out to a certain maximum 
i which increases with redshift due to projection effects 
(see Fig.EJ. The effect of losing the lowest multipoles due 
to survey boundaries or systematic effects or the highest 
redshifts due to depth and photometric redshift errors is 
thus more substantial for distinguishing between cluster- 
ing models (see Fig.lHJ. 

These results remain robust when varying the fiducial 
model in w if interpreted as the significance of detecting 
a given variation in the gravitational potential due to 
dark energy clustering on scales corresponding to c^r\Q ~ 
14cg Gpc. We show this rescaling of Fig. El to S/N ~ 1 
for effective la errors on the change in the gravitational 
potential due to dark energy clustering out to a given 
scale. The caveat is that as — 1 one would predict 
no change out through the horizon in an adiabatic model. 
These results are also robust to changes in the halo model 
for galaxy bias, variations in the galaxy source density, 
and improvements in photometric redshift accuracy. 
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V. DISCUSSION 

The ISW effect provides a unique probe of the dynam- 
ics of the dark energy through its high sensitivity to the 
smoothness of the dark energy. Although the significance 
of its detection through cross correlation with galaxy sur- 
veys or gravitational lensing statistics will never exceed 
^ 10a due to the cosmic variance of the CMB tempera- 
ture field , it is likely to be the best means of probing 
the clustering dynamics since their effects appear only on 
the largest scales. 

A 10(7 detection of the correlation will allow its mea- 
surement to ^ 10% and hence yield the ability to limit 

3% variations in the gravitational potential due to the 
dark energy clustering if they appear on the Gpc scale. 
In the adiabatic model, this would allow for a test of the 
quintessence, or canonical single scalar field, hypothesis 
if |l + w| ^ 0.05. For distinguishing clustering models, as 
opposed to a simple detection of the correlation, it is im- 
portant to measure correlation at the largest angles out 
to z ~ 2. That will require a deep nearly all-sky survey 
such as LSST [2^. Furthermore the galaxy clustering 
signal on very large angles is itself small and controlling 



systematics in the galaxy surveys will pose a significant 
challenge (e.g. pll27t). 

Finally, beyond the adiabatic model, clustering in the 
dark energy may have strong effects even ii w — > — 1 . For 
example, in the isocurvature models designed to suppress 
the quadrupole [2ll , the change in the gravitational 
potential is about 6 times that of the adiabatic model. 
Such models are potentially testable with CMB-galaxy 
cross correlations in the right redshift range. Likewise 
explanations for the acceleration that involve modifica- 
tions to gravity in place of the dark energy may predict 
a different CMB-galaxy correlation on the largest scales. 
The uniqueness of this probe of the acceleration of the ex- 
pansion should motivate future studies despite challenges 
facing its exploitation. 
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